Epigenetics & Chromatin

Schistosoma mansoni is a parasitic flatworm that has two, genetically determined, sexes. We used aggregated data of 8 posttranslational histone modifications (ChIP-Seq), chromatin accessibility (ATAC-Seq), transcription (RNA-Seq) and genome feature annotations to decipher the histone code of genes involved in the differentiation of schistosome gonads (i.e. female ovaries and male testes). We show schistosome gonads express at least two classes of protein coding genes: H3K4me3-positive genes that display canonical features of eukaryotic protein-coding genes such as peaks of H3K4me3 at the transcription start sites (TSS) and increases in histone acetylation marks towards the transcription end site (TES), but also a non-canonical H3K9/27me3 plateau just upstream of the TSS. H3K4me3 enrichment at the TSS is highly predictive for transcription strength in these genes compared to a second class of protein coding genes (H3K4me3-negative genes) that do not display this pattern and is characterised by absence of the investigated histone marks at TSS and TES. This is indicative of the existence of hitherto unknown, potentially schistosome-specific histone marks in these genes. The absence of H3K4me3 at the TSS is not associated with inducible or stable gene expression in the gonads. Instead, gene ontology analysis indicates that H3K4me3-positive genes are related to functions which typically govern processes such as metabolism or signal transduction while H3K4me3-negative genes are dedicated to cell communication or immune responses. Second, individual histone modifications and their combinations are associated with functional features of the schistosome genome, known as "chromatin colours". In H3K4me3-positive genes, there is clear co-linearity of 3 colours, which strongly suggests a functional role for histone modifications in the control of transcription pre-initiation, promotor release, and transcription termination. Third, there are striking chromatin structure changes during maturation of the gonads in all genomic features including protein-coding and non-protein coding genes as well as repetitive sequences. The nature of these changes is different in both sexes. H3K36me3 and H3K9me3, as well as H3K23ac and H3K9ac show the strongest variations. Last, we show that pharmacological inhibition of histone demethylation activity by IOX1 leads to a concentration-dependent separation ("divorce") of schistosome couples confirming the importance of H3K36/H3K9 methylation for pairing maintenance and indicating histone demethylases as a potential drug target family. Collectively, our findings offer unprecedented insights into histone codes and chromatin dynamics governing the reproductive development of S. mansoni gonads.

BACKGROUNDDuring spermiogenesis, histones are exchanged by protamines (PRMs) in spermatids, which results in DNA hypercondensation and protection. Rodents and primates express two PRMs (PRM1 and PRM2) in a species-specific ratio. Maintaining this ratio is necessary for functional chromatin reorganization and alteration is associated with sub- or infertility in mice and humans. Prm1 and Prm2 deficient mice are infertile, while Prm1+/- males are subfertile showing a severely altered PRM ratio. Prm2+/- males are fertile and display a protamine ratio comparable to WT. OBJECTIVESHere, we addressed the question whether loss of one allele of Prm1 and one allele of Prm2 affects fertility. MATERIAL AND METHODSDouble heterozygous (dHET) mice lacking one allele of Prm1 and one allele of Prm2 were generated and analyzed RESULTSdHET males were infertile with sperm showing retention of histones and TNPs, high levels of PRM2 precursor and decreased levels of mature PRM2. In mature sperm the PRM ratio and the total PRM content was not altered. However, CMA3 staining revealed incomplete protamination and sperm nuclei appeared more rounded and slightly bigger, suggesting impaired DNA-hypercondensation. In dHET sperm, DNA degradation was apparent, but to a lower level compared to sperm from Prm1 and Prm2 deficient males. Increased 8-OHdG levels suggested oxidative stress in the epididymis of dHET mice. However, a fraction of dHET sperm were capable of fertilization, with embryonic development up to 8-cell stage. DISCUSSION AND CONCLUSIONThese results suggest, that male factor infertility might not be reliably detected by measuring PRM1/PRM2 ratio but rather by determining the level of protamination by e.g. CMA3 analysis and pre-PRM2 retention.

The SWI/SNF complex comprising of the catalytic subunit, Snf2, is a key regulator of gene expression and DNA damage repair in eukaryotic cell. Candida albicans Snf2 is known to regulate hyphal formation. In this paper, we have investigated the role of this protein in DNA damage response. We show that CaSnf2 is required for cell division as deletion of both copies of SNF2 leads to increased duplication time. The mutant cells form clumps with increased chitin and {beta}-glucan deposition on the cell wall. The altered cell wall phenotype leads to reduced uptake of genotoxic stressors leading to increased resistance to both methyl methane sulfonate (MMS) and hydroxyurea (HU). In addition, resistance of Casnf2{Delta} cells to MMS also appears to be mediated by upregulation of CaRAD9 expression by CaFun30, an ATP-dependent chromatin remodeling protein, and CaRtt109, a fungal-specific histone acetyltransferase. The response of Casnf2{Delta} to genotoxic stressors is at variance with the response of Scsnf2{Delta} mutant, highlighting the differences in DNA damage response/repair pathway between the two organisms. Finally, we show that Casnf2{Delta} mutants are extremely sensitive to azoles due to downregulation of multi-drug resistance pumps leading to reduced efflux of the drug.

BackgroundDDX53 (DEAD-box helicase 53, known also as CAGE) is an intronless gene on the X chromosome, which expression shows strong testis specificity. It belongs to the group of cancer-testis (CT) antigens, with most studies to date focusing on its role in cancer, but the precise biological function of DDX53 remains unclear. Previous reports identifying rare DDX53 variants in infertile men provided the rationale for investigating the role of DDX53 in the context of human spermatogenesis. By using the human seminoma cell line (TCam-2) as an in vitro male germline model, we aimed to investigate the function and molecular targets of DDX53. MethodsIn our study, we used transcriptomic and proteomic approaches (RNA sequencing (RNA-seq), enhanced crosslinking and immunoprecipitation (eCLIP), and Co-immunoprecipitation coupled with Mass Spectrometry (Co-IP-MS)) to investigate the role of DDX53 in the context of human spermatogenesis. By using modified TCam-2 cells to express either DDX53-FLAG or GFP-FLAG, we identified regulated genes, RNA targets, and potential protein interactors of DDX53. In addition, we employed Western Blot, RT-qPCR, immunostaining, and confocal microscopy to gain deeper insight into the DDX53 protein. ResultsOur RNA-seq and eCLIP data provide evidence that DDX53 regulates gene expression changes and directly interacts with a broad spectrum of RNA transcripts. Moreover, for the first time, we described RNAs and protein interactors of DDX53 in the context of spermatogenesis. Subcellular localization analysis by confocal microscopy indicated a predominantly cytoplasmic distribution of DDX53, with partial nuclear presence in TCam-2 cells. We also identified DDX53-positive structures that may correspond to germ granule-like assemblies, although their precise nature remains to be determined. Additionally, we confirmed DDX53 presence in human testis using a specific, commercially available anti-DDX53 antibody. ConclusionsThis studys data indicate that DDX53 protein acts as a regulator of RNA metabolism in human cells. Collectively, we show that it participates in transcriptome regulation (including splicing) in male germ cells and exhibits transcriptome-wide RNA interactions, but its wider biological role remains to be clarified.

Nuclear blebs are herniations of the nucleus that occur in many human conditions including aging, heart disease, muscular dystrophy, and many cancers. Nuclear blebbing causes nuclear rupture and cellular dysfunction. However, understanding the formation, stability, and identification of nuclear blebs remains an ongoing challenge. Our previous studies reveal that nuclear blebs are best hallmarked by decreased DNA density. To determine if chromatin decompaction underlies decreased DNA density in nuclear blebs, we investigated the histone composition of nuclear blebs across multiple cell lines. Time lapse and immunofluorescence imaging revealed that global histone H2B and H3 levels are decreased in the nuclear bleb relative to the nuclear body. Next, we imaged histone modification states of euchromatin and heterochromatin, which respectively track decompact and compact states of chromatin. Overall, we find that nuclear blebs display variable histone modification state across cell lines, as euchromatin does not consistently enrich nor is heterochromatin consistently depleted. Nuclear blebs did consistently show active RNA Pol II initiation is enriched relative to elongation. Thus, we find that the local histone modification state is not an essential component of nuclear blebs while transcription initiation enrichment over elongation is reproducible across cell lines and conditions. Summary statementWe measured histones and their modification states in nuclear blebs. We find that chromatin state is variable while transcription initiation is consistently enriched relative to elongation in nuclear blebs.

4C-seq is a cost-effective 3C-based assay that measures the interactions between a single genomic element and all other genomic elements. However, 4C-seq data remains semi-quantitative because it cannot be deduplicated without UMIs. To address this, we developed an open source method, FourC, based on a Bayesian Bernoulli regression model, that overcomes the duplication problem and models spatial patterns with Gaussian processes to identify significantly enriched and differential contacts. We demonstrate the utility of FourC on 4C-seq data that profiles the local chromatin structure at key genes necessary for pancreatic differentiation and under CRISPR perturbation of enhancers.

Epimutations are changes in chromatin modifications, such as DNA methylation or histone modifications. Some of these epigenetic changes can be inherited for several generations, and they potentially contribute to evolutionary processes. Estimates of epimutation rates now exists in a few species, but the presence and function of epigenetic marks are not conserved across different species. To understand the properties of epimutations in fungi, we performed a mutation accumulation experiment with the filamentous fungus Neurospora crassa and investigated spontaneous changes in DNA methylation and trimethylation of lysine 9 on histone H3 (H3K9me3) in the mutation accumulation lines. We observed that centromeric regions are hotspots of spontaneous DNA methylation changes in N. crassa. In these hotspot regions, DNA methylation changes were transmitted across mitoses, but changes occurring in euchromatin were not maintained. The rate of DNA methylation changes was around 30 000 fold faster than the genetic mutation rate. We did not observe spontaneous changes in H3K9me3 that were transmitted across mitoses. Our results show that while spontaneous epimutations occur in this species, they occur predominantly in gene poor heterochromatic regions, so their impact for evolutionary adaptation may be limited.

Purine moieties are essential for many functions within the eukaryotic cell, including energy, signaling and nucleic acid synthesis. While purine starvation is known to induce stress resistance in eukaryotic model organism budding yeast Saccharomyces cerevisiae, it remains unclear whether the physiological response is related to disruption of synthesis pathway in particular position or it is uniform across all genetic deficiencies within the de novo adenine biosynthesis pathway. It is also not known how purine starved cells perceive purine shortage - weather they share the same signaling elements with nitrogen starvation or not. MethodsWe characterised physiology of strains with deletions in adenine biosynthesis pathway when cultivated in full or purine deficient and compared to cell physiological parameters when cultivated in nitrogen deficient media. We tested stress tolerance, carbon flux, cell cycle arrest and did transcription profiling (RNA-seq). ResultsOur findings demonstrate that purine starvation-induced stress resistance is significantly modulated by the specific step at which the pathway is interrupted. Transcriptional analysis revealed that purine starvation in many aspects phenocopies nitrogen starvation, particularly - in both starvations strong downregulation of ribosome related genes occurs. In the same time several metabolic features which differ from N- and ade- starvations: pentose phosphate pathway is specifically upregulated within ade4{Delta}-ade2{Delta} and downregulated in N-cells. Notably, the expression of stress-responsive genes such as HSP12, HSP26, and GRE1 varied between mutants, suggesting that the accumulation of pathway intermediates (e.g., AIR in ade2{Delta}) or the absence of downstream precursors (AICAR) alters the perception of starvation especially in the case of ade16{Delta}ade17{Delta} strain. ConclusionsMetabolic and stress-tolerance phenotypes of purine auxotrophs are not merely a result of purine depletion but seems that the response is signalled via the same pathways, like TOR1. The results suggest that strains having mutations within various positions of the purine pathway "perceive" purine limitation a bit differently - especially when we compare the end of the pathway with the other mutants. Different phenotypic outcomes of the occasional purine depletion might give preferences for organisms which have mutations in the beginning rather at the end of the pathway. Besides, our findings might have implications in the design of synthetic pathways and the use of auxotrophic markers in yeast research.

DNA methylation is a key epigenetic mechanism influencing gene regulation and cellular identity. In skeletal muscle, methylation contributes to fiber-type specification, metabolic programming, and satellite cell function, with evidence of sex-specific differences. Here, we investigated whether spatial regionalization of gene expression along the proximal-distal axis of the tibialis anterior (TA) is mirrored by corresponding patterns of DNA methylation. Using MeDseq on TA sections from muscles previously analyzed by spatial transcriptomics, we profiled methylation across transcriptional start sites (TSS), gene bodies, and regulatory elements. Despite robust spatial differences in transcriptomes, methylation patterns were largely uniform along the proximal-distal axis, indicating that DNA methylation does not underlie regional gene expression in adult TA muscle. In contrast, sex emerged as the primary determinant of methylation variation. Male muscles exhibited widespread hypermethylation at TSS, gene-bodies and regulatory regions, corresponding with sex-specific transcriptional programs, including glycolytic fiber enrichment in males and oxidative fiber markers in females. Notably, chromatin- and methylation-associated regulators such as Setd7, Gsk3a, and Bmyc were upregulated in males, suggesting mechanisms linking transcriptional control to epigenetic state. These findings highlight that while spatial gene expression is transcriptionally driven, sex-specific epigenetic programs dominate adult skeletal muscle, underscoring the need to consider sex in multi-omic studies of muscle biology.

Epigenetic mechanisms regulate gene-expression by altering the structure of the chromatin without modifying the underlying DNA sequence. Histone post-translational modifications (PTMs) are critical epigenetic signals that influence transcriptional activity, promoting or repressing gene-expression.Understanding the impact of individual PTMs and the combinatorial effects is essential to deciphering gene regulatory mechanisms.In this study,we analyzed the ChIP-seq data for 26 PTMs in yeast, examining the PTM intensities gene-wise from positions-3 to 8 in each gene.Using XGBoost classifiers, we predicted gene transcription rates and identified key histone modifications and nucleosomal positions that are critical in gene-expression using explainability measures (such as SHAP). Our study provides a comprehensive insight into the histone modifications, their positions and their combinations that are most critical in gene regulation in yeast.The proposed explainable Machine Learning models can be easily extended to other model organisms to provide meaningful insights into gene regulation by epigenetic mechanisms.

BackgroundThe histone methyltransferase EZH2, enzymatic core of the trimeric polycomb repressive complex 2 (PRC2), has been shown to promote small cell lung cancer (SCLC) survival through epigenetic silencing of multiple targets including Class I MHC molecules (HLA-A/B) and DNA repair factors (SLFN11). Treatment of SCLC cells with EZH2 inhibitors in vitro can reactivate expression of these genes and result in therapeutic response to immune checkpoint inhibition (ICI) and chemotherapy. Here, we investigate the impact of EZH1/2 dual inhibition on 3D chromatin structure and its relationship to transcriptional regulation in neuroendocrine (NE) SCLC. ResultsEmploying Micro-C, a micrococcal nuclease-based 3D genome mapping technique, we show that EZH1/2 inhibition with Valemetostat induced significant changes at multiple genome organizational levels (compartment, topological associated domain, and chromatin loop) without incurring cell death in NE SCLC. Alterations in 3D genome permissive for transcriptional activation were correlated with increased chromatin accessibility (ATAC-sequencing) and expression of target genes (transcriptome profiling). Known transcription factor motif discovery revealed enrichment of non-NE motifs (e.g., REST) in regions with gained chromatin accessibility in Valemetostat-treated cells, consistent with results from gene set enrichment analysis demonstrating NE to non-neuroendocrine lineage shift. Notably, EZH1/2 inhibition reactivated Class I MHC expression by facilitating enhancer-promoter looping. ConclusionOur results demonstrate that repression of a subset of EZH2 targets including Class I MHC genes is affected through modulation of 3D genome structure to the level of chromatin looping and further support clinical investigation of EZH2 inhibition in boosting therapeutic efficacy of ICI in SCLC patients.

BACKGROUNDDNA methylation is the most prevalent epigenetic modification in human cells and undergoes dynamic changes during cell differentiation, disease progression, and aging. Here, we introduce Locus-Enriched Mapping Of Nucleotide methylation (LEMONmethyl-seq): an optimized, cost-effective pipeline for single-nucleotide detection of DNA methylation using locus-specific amplification and long-read DNA sequencing. RESULTSWe apply LEMONmethyl-seq to profile DNA methylation of endogenous gene promoters across different cell types along with DNA methylation establishment and long-range propagation induced by CRISPR epigenome editing technologies. We profile dynamic changes in DNA methylation patterns on transposable element genomic loci during global epigenetic resetting in stem cells, and we identify site-specific enrichment of non-canonical CpH methylation on genomic sites in stem cells and cultured neurons. Lastly, we apply LEMONmethyl-seq to profile DNA methylation across the MGMT promoter, a clinical biomarker for glioblastoma. We identify additional differentially methylated sites correlated with chemotherapeutic sensitivity, which may be clinically relevant. CONCLUSIONSTogether, LEMONmethyl-seq serves as a cost-effective, long-read DNA methylation sequencing pipeline that advances methods for detecting DNA methylation patterns and dynamics in mammalian cells. We envision its broad use for studying chromatin pathways, diagnostics, and therapeutic applications.

Escherichia coli is highly sensitive to acid and osmotic stress but adapts by modulating the expression of stress responsive genes. Nucleoid-associated proteins (NAPs) play key roles in DNA organization and sensing environmental changes. The histone-like nucleoid structuring protein H-NS is an NAP acting as a global regulator of stress genes. H-NS may alter local chromatin structure to modulate the expression of such genes in response to environmental stress. The H-NS homolog StpA co-regulates several target genes, but its precise role is poorly defined. To investigate the regulatory interplay between these two proteins, we examined transcription, DNA binding and chromatin structure at two regulated operons, hdeAB and proVWX, in E. coli following exposure to acid and salt shock. Our results show that H-NS senses pH and osmotic cues to remodel chromatin and relieve repression, while StpA compensates for H-NS loss, particularly at proVWX, highlighting a coordinated regulatory network.

During early embryonic development, cells transition from naive to primed pluripotent state. Various culture conditions have been established to revert primed cells back to naive state, to increase differentiation potential and to reset epigenetic abnormalities. In this study, we modified culture conditions to allow primed-to-naive conversion under feeder-independent and normoxic conditions (FINO medium), which exemplified the need for a quantitative measure of pluripotent states. DNA methylation (DNAm) profiling revealed extensive hypomethylation at naive state, but also significant gains of methylation at specific sites in the genome. We demonstrate that DNAm patterns can be used to benchmark culture protocols. Furthermore, we developed a naive-score based on DNAm at two genomic sites, which can be analyzed by digital PCR to monitor transition between pluripotent states. Our study describes a simplified culture protocol for primed-to-naive conversion, offers insights into the specific DNAm changes, and introduces a robust DNAm-based biomarker to track this process effectively.

BackgroundMaize knobs are regions of constitutive heterochromatin that are readily identified in both meiotic and somatic chromosomes. These structures have been characterized as stable throughout the cell cycle, exhibiting late replication during the S-phase, and are composed of two specific families of highly repetitive DNA sequences: K180 and TR-1. Although widely used as cytogenetic markers due to their variability in number and chromosomal position across inbred lines, hybrids, and landraces, little is known about their chromatin structure and dynamics. In this study, we analyzed chromatin accessibility of knobs using DNS-seq data across four maize tissues representing distinct developmental stages. ResultsOur results reveal that K180 knobs exhibit tissue-specific variation in chromatin accessibility, transitioning between open and closed states during development. In contrast, the TR-1 knob of chromosome 4 remained consistently inaccessible across all tissues analyzed. A knob composed of both K180, and TR-1 further supported this observation, with only the K180 region showing dynamic accessibility. To validate these findings, we also analyzed other repetitive regions such as centromeres, which showed a uniformly closed chromatin structure similar to TR-1. These results suggest a unique developmental modulation of chromatin accessibility associated with K180 repeats. While the chromatin accessibility of knobs does not reach the levels observed at Transcription Start Sites (TSS), the comparison among different classes of repetitive DNA within maize constitutive heterochromatin provides compelling evidence for sequence-specific and tissue-specific chromatin dynamics. ConclusionsOur findings uncover a previously unrecognized property of maize knobs and establish a reference for future studies on chromatin organization and epigenetic regulation of repetitive DNA in plant genomes.

The human genome is partitioned at different levels of 3D genome organization, with topologically associating domains (TADs) being among the most well-known and biologically important structures. TAD boundary disruption is associated with a wide range of diseases such as cancer, neurological and developmental disorders. Numerous methods have been developed to detect TAD boundaries from chromatin contact maps obtained with Hi-C technology. However, these methods are largely limited by the resolution of Hi-C data, typically 1 Kb to 100 Kb. In contrast, functional DNA loci, collectively referred to as epigenomic data, are profiled at a much higher resolution (100-200 bp for a typical ChIP-seq experiment). To improve the resolution of boundary detection, we hypothesize that the patterns of epigenomic signals associated with regions in proximity to TAD boundaries can serve as embeddings for these genomic regions, defining region similarity. These embeddings, along with their positional relationships, can be effectively modeled using deep learning to achieve more precise boundary prediction. We present EpiTADformer, a transformer-based model that takes as input transcriptional and histone modification signals of neighboring regions centered around TAD boundaries. We demonstrate that EpiTADformer outperforms feedforward neural network, convolutional neural network (CNN), and bidirectional long short-term memory (BiLSTM) network architectures. These results suggest the positional information of epigenomic signals surrounding TAD boundaries provides a strong predictive signal, enabling improved performance of the transformer model. Our findings highlight the potential of epigenomic signals to serve as region embeddings for refining the epigenomic language of TAD domains and 3D genome organization.

Ehmt2 is a key H3K9 methyltransferase that regulates genome silencing and structural integrity during animal development. In addition to this canonical function, Ehmt2 has also been implicated in neural tissues mediating both direct and indirect transcriptional activation, and exon splicing, to facilitate proper neural cell differentiation and survival. Several germline loss-of-function animal models have been developed showing both conserved and divergent phenotypes that range from embryonic lethality to behavioural deficits in adult, fertile animals. Here, we generated the first maternal-zygotic ehmt2 loss of function mutant in zebrafish using CRISPR-Cas9 mutagenesis. An assessment of the pattern of H3K9 methylation in mutant embryos by ChIP-seq indicates that there are aberrant levels of this repressive mark, including reduction in discrete 5 non-coding regions of genes, but with no significant change in the overall pattern distribution of these marks across the genome. Global transcriptome and morphological analyses demonstrated that mutant embryos displayed greater variation in the timing of developmental progression that is, on average, slower compared to controls. Despite this, mutant embryos ultimately survive and are fertile. Through examination of progenitor cell dynamics and gene expression profiles, we found that the delay in embryonic development was associated with longer rates of S-M phases of the progenitor cell cycle in mutants leading to deficits in tissue growth. Finally, our data suggest a robust network of epigenetic regulators can potentially compensate for Ehmt2 loss of function and permit embryonic development and survival in ehmt2 mutant zebrafish. Our work establishes a zebrafish ehmt2 loss of function model that will facilitate examination of the complex and varied roles of Ehmt2 in vertebrate development.

Primordial germ cells (PGCs) are the population of cells that, in the human embryo, specify day 12 post-fertilization, and form the precursor cells for the future egg or sperm cells. Although in vitro differentiation of PGCs from human stem cells has been achieved, these primordial germ cell-like cells (hPGCLCs) fail to further mature. The reason for this is unclear. Previous studies in mice revealed that several specific metabolic changes occur during the maturation of these cells, which are essential for their developmental progress. However, very little is known about the metabolic profile of human primordial germ cells. In the severe scarcity of human PGCs, hPGCLCs serve as a research model to study PGC formation. To investigate this, we differentiated hPGCLCs using induced-pluripotent stem cells and performed a mass spectrometry analysis to establish their metabolome and proteome. These cells revealed distinct metabolic profile, with changes particularly at the proteome level. This included a shift between canonical and non-canonical citric acid cycle in hPGCLC, downregulation of late-stage glycolysis and reduction of nucleotide de novo synthesis. By providing an integrative map of these metabolic networks, we aim to provide insight on the influence of metabolism on human PGC development that could help improve methods for in vitro differentiation and maturation hPGCLCs.

Sirtuin 6 (SIRT6) is an important regulator of DNA repair, metabolism, chromatin maintenance and longevity. SIRT6 Serine 10 phosphorylation controls SIRT6 recruitment to the sites of DNA damage. To explore the effect of SIRT6 Serine 10 phosphorylation on lifespan, we generated two SIRT6 mutant mouse strains: phospho-null S10A and phosphomimetic S10E. The S10E mutant mice demonstrated enhanced DNA repair capacity, elevated LINE1 expression and reduced lifespan in male mice compared to the wild-type and S10A mice. This result suggests that SIRT6 S10E mutation enhances DNA repair capacity at the expense of reduced LINE1 silencing leading to shorter lifespan. While both SIRT6 functions in DNA repair and chromatin maintenance are important for longevity, our results suggest that when the balance between these functions is shifted, diminished of LINE1 control has a stronger impact on lifespan than enhanced DNA repair.

DNA methylation is a regulator of bacterial gene expression and adaptation, influencing traits such as virulence and antimicrobial resistance. The dynamic nature of DNA methylation enables rapid responses to changing environments and is a source of heterogeneity in bacterial populations. However, condition-dependent DNA methylation and consequences for transcriptional output remain poorly understood. We applied Oxford Nanopore sequencing to profile DNA methylation during exponential growth and late stationary phase of Salmonella enterica serovar Typhimurium and integrated these data with transcriptomic analyses. We found that each DNA methyltransferase (MTases) exhibits a distinct activity pattern across growth stages, which could not be explained by transcriptional levels of the corresponding enzymes. As predicted, DNA methylation patterns determined by regulatory MTases were dynamic across growth conditions whereas methylation patterns of MTases belonging to R-M systems were comparatively stable. We identified growth stage-specific methylation patterns for all studied MTases and correlations between methylation states and gene expression patterns. Together, these findings chart DNA methylation networks in the epigenetic regulation of bacterial physiology. Author summaryDNA methylation in bacteria is best known for its role protecting DNA from endonucleases, such as restriction-modification, and coordinating chromosome replication and mutation repair, yet DNA methylation also regulates gene expression and cell physiology. Previous studies primarily examined bacterial DNA methylation at single time points or in limited genomic regions, providing only a partial view of its biological significance. In this study, we used Oxford Nanopore sequencing to compare DNA methylation patterns in Salmonella enterica during exponential growth and late stationary phase then integrated these data with corresponding gene expression profiles. We identified numerous methylation target motifs, all of which demonstrated constitutively methylated or unmethylated regions. This systems-level analysis clarifies the role of DNA methylation in bacterial adaptation across growth stages and demonstrates the utility of Oxford Nanopore sequencing for genome-wide methylation profiling.